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Abstract 

The aim of this lecture is to present in a comprehensible way what the 
BRST quantization means and how the "classical" master equation, action 
and BRST transformations have to be prolonged towards the same "quan- 
tum" items. The presentation will focus not only on the standard BRST 
symmetry, but on larger symmetries as sp{2), both in the Lagrangean and 
in the Hamiltonian formalisms. How to find answers to these questions in 
more sophisticated cases will be illustrated by the example of a nonlinear 
system with open superalgebra. 

1 Introduction 

All fundamental interactions from nature are described by theories with internal 
symmetries or gauge theories. Such theories, as QED, QCD, electro- weak theory 
and gravity, ask for a special approach because of the unphysical degrees of 
freedom involved by the gauge invariance. For example, when quantizing such 
theories the direct computation of the path-integrals is meaningless since the 
integration over gauge directions in the measure would make the path-integral 
infinite- valued. The BRST technique [T] is one of the approaches which allows 
to overcome this dificulty. It assumes the replacement of the local (gauge) 
symmetry with a global (BRST) one. The BRST symmetry is expressed either 
as a differential operator s, or in a canonical form, by the antibracket ( , ) in the 
Lagrangean (Batalin-Vilkovisky) case |2] and by the extended Poisson bracket 
{ , } in the Hamiltonian (Batalin-Fradkin-Vilkovisky) formulation [3j: 

s* = (*,5) = {*,f]}. (1) 

The BRST generator S and the BRST charge O are both defined in extended 
spaces generated by the real and by the ghost-type variables. To perform path- 
integral calculations in this frame, it is necessary to remove the redundant 
gauge variables, that is to gauge-fix the action by choosing a gauge fermion 
^.Elimination of gauge variables assures the BRST invariance of the action [3] 
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but not of the measure. A BRST transformation of the coordinates could gen- 
erate non-trivial terms in the action. These terms can be exponentiated and 
adsorbed in the action ( Fadeev- Popov trick). One obtain an extended action 
called "quantum" action. 

In this paper we will explicitly show how the BRST formalism allows to con- 
struct well defined path-integrals and what is the concrete form of the quantum 
master equation in the special case of a physical system described by a nonlin- 
ear gauge algebra. The interest on such systems was initiated by the discovery 
of conformal field theories [S] which led to a new class of gauge theories with 
nonlinear gauge algebras, the so-called Wat algebras [5]. The BRST approach 
for non-linear superalgebras was developed in [7] . 

In the standard BRST approach the gauge fixing procedure asks for a " non- 
minimal sector" which suposes the introduction of new, supplementary vari- 
ables. In the extended symmetries this non-minimal sector appears in a natural 
way. Now, the gauge fixing function can be constructed in a more simple way 
and the gauge fixing action is obtained easily. Singura difficulty which can be 
appear in the case of the extended BRST formalism is asigurarea of the invari- 
ance of the integrating measure, so that the average of the observables to be 
independent de the choice of the gauge. Ne propunem in this paper to show how 
can done this in the case of gauge theories with nonlinear gauge superalgebras. 



There arc two main procedures in obtaining a quantum description of the gauge 
systems: the Dirac quantization method and the path integral method. In the 
Hamiltonian formalism the existence of the gauge symmetries is equivalent with 
the existence of some constraints imposed to the theory and this is why the 
Dirac quantization method, which is especially apropriate for the description of 
systems with constraints, is a very useful method in quantizing gauge theories. 
As we intend to built up the quantum equation for interacting fields, we will 
concentrate our approach on the path integrals' formalism which is more suitable 
for the quantum description of these systems. This approach to the quantization 
of systems with symmetries, both in the Lagrangean (BV) formalism and in the 
Hamiltonian (BFV) one, started from the Fadeev-Popov trick and ended with 
the development of the BRST technique. The gauge fixing procedure assures 
the BRST invariance of the gauge fixed action but not automatically of the 
generating functional Zy and of the expectation values {F)y- For a gauge 
theory with an action S[(l), </>*] we have: 



2 Quantum Master Equations 



2.1 The Standard BRST Approach 





and 




(3) 
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The independence on the choice of Y and the BRST invariance of ^ and ^ , 
essential requirements for a consistent quantum description, asks more for the 
invariance of the measure d<j)^ to the BRST transformation: 

^ ^^A^ s)^i - 0^ + (4) 

o<Pa 

where /i is a constant, anti-commuting parameter. In the general case the 
transformation (j4]) leads to 

= (1 - ^ (1 - (A5)M)d0^ (5) 

5(1) <5<Pb 

and nontrivial terms in the measure could be generated. The invariance condi- 
tion would impose: 

d0'^ = d0^^A5-O;A^(-ir--^-|^. (6) 

The condition ([6]) might not be fulfilled. Although, all the terms generated by 
(AS')/i)(i0"^ can be exponentiated and absorbed as "quantum corrections" of 
the action S: 

W + nWi + TiVa + • • • (7) 

One obtains the " quantum action" W which satisfies a Quantum Master Equa- 
tion (QME), an extension of the classical master equation (S", S) — 0. It has the 
form: 

i(I^,VK) -i^iAPF = 0. (8) 
The equation ([5]) can be written in the equivalent form 

Ae^*^ = 0. (9) 

It is important to note that A is not a true derivation, its action beeing given 
by the following rule: 

A(a/?) = aA/3 - (~)""(Aa)/3 - i?iA(a,/3). (10) 

The classical BRST operator s can be, it also, extended to a quantum one: 

s* = (*, S) ^ a* = W) - ih* (11) 

The quantum version of the BRST invariance will be expressed by: 

aF = ^ A(i^e^'^) = 0. (12) 
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2.2 The sp{2) extended BRST Quantization 

Because of some difficulties the standard BRST approach met in the gauge fix- 
ing procedure, procedure which suposes the introduction of some supplementary 
variables from a non-minimal sector, extended BRST symmetries has been for- 
mulated. The main such extension is known as the sp(2) BRST symmetry and 
suposes the existence of two anticommuting differential operators, si and S2, 
which can be joined in a symplectic doublet si,S2 with: 

s = si -t-S2; s^ = 0. (13) 

In the Hamiltonian formalism [S], the sp{2) BRST symmetry is canonical 
generated by the BRST charges ila, a = 1,2 with e{il.a) = 1- By defining an 
extended phase-space and a generalized Poisson bracket, the BRST charges will 
be given by: 

Sa* = {*,na};a = 1,2. (14) 

Their concrete form depends on the theory and it can be obtained by using 
the homological perturbations theory [4]. The extended Hamiltonian H is the 
solution of the equations 

{H,n"'} = 0, a = 1,2 

with boundary condition 

H\q^x=o = Ho{q,p). 

It is obtained using the same homological perturbations theory 4 . The gauge 
fixing procedure leads to a gauge fixed Hamiltonian of the form 

where Y is the gauge fixing functional defined in terms of the real and ghost 
type coordinates. 

Even if it is easier to develop the Hamiltonian formalism, it does not always 
lead to a covariant gauge fixing action. This is way the passage to the La- 
grangean formalism is indicated. It is in the Lagrangean frame where we should 
study the invariance of the generating functional Zy and of the expectation 
values of the form ([3]) for observables. 

In the Lagrangean case [9], two antibrackets are defined and the master 
equations take the form: 

^iS,S)a + VaS = 0; a = 1,2 (15) 

where {Va,a = 1,2} represent two special "non-canonical" operators. There 
is a direct modality of obtaining the sp{2) BRST Lagrangean formalism, but 
it assumes the use of a very large spectrum of ghost generators. To avoid this 
unuseful extension, it is simpler at the classical level to construct the Lagrangean 
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formalism following its equivalence with the Hamiltonian one [TU]. We will 
present how the classical sp(2) theory can be extended at quantum level in 
the frame of path integrals formalism. We will not give here more details but 
we direct to where quantum sp(2) master equations were obtained in the 
Hamiltonian and then in the Lagrangean schemes. In this quantum context, the 
invariance condition ^ has as correspondent the equations: 



A'^S = 0; A" = (-1)^-^ — -^ + v'' = (-1)^-^ — -^ + e''''^Ab^ (16) 



and 

^a^b ^ ^b^a ^ ^ ^ -^ 2. (17) 

If A^S* ^ 0, the operators A*", a = 1, 2 determine the sp{2)- QME in the sp{2) 
symmetric formulation of gauge theories: 

A^ei^ = {)^'^{W,WY = inA''W, a = 1,2 (18) 

where W represents the "quantum action", the ?i— order corrections being in- 
volved by integrating measure. In the next section we will come back to this 
problem, approaching it in more concrete and applied manner. 



3 The Example of Nonlinear Superalgebras 

Let us consider a system described in a phase-space M = {q^,pi,i — 1, • • • ,n} 
with e{q^) ~ e{pi) = by a set of first class constraints {Gq, a = 1, • • • , m} 
with s{Ga) — Sa and by the canonical action 

Scan[q,P,u]^ j dt[qpi- H^°\q,p,u)], i = 1,- ■ ■ ,n (19) 

= Ho(g,p)+""G„. (20) 

The Lagrange multipliers {u",a = I,-- - ,m} will play the key rol in the es- 
tablishment of the equivalence between the extended BRST Hamiltonian and 
Lagrangean formalisms. We suppose that the constraints satisfy involution re- 
lations of the form: 

{G^,Gp}^G^flp + G,Gsgl'p (21) 
The gauge algebra is given by the relations (|2T1) and by 

{i?o, G„} = GpVi + GsGpU!^'. (22) 

We consider the case when f2p, 5^^, , U^^ are true constants. We deal in 
this case with a quadratic non-linear Lie algebra, case which correspond to a 
constrained system with open algebra [12] ■ An example of such system is the 
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Yang-Mills fields theory with the nilpotent BRST charge a quadratic function 
of the Fadeev-Popov ghost fields. 

We will develop the sp(2) BRST Hamiltonian formalism for the previous 
nonlinear superalgebra. The extended phase space will be generated in this 
case by the real variables {q^,Pi,i = l,...,n} and by the ghost- type variables 
{Q"", 7^/36, A", TT^}. The last ones satisfythe relations: 

On the basis of the Jacobi identities for the structure functions of the superal- 
gebra we can show that for any superalgebras (PT|) . ([^ t 

Hipp + Gsg%){f2c. + Gsgil){^)'^''+'^^'^+'-^+'-'^ + 

+ (/^/j + Gsg%){f2p + G,-.9^p(~)(^^+i)(^°+^^)+^^^i£,,^^Q^'^Q'"'Q-. (23) 

If in (jl9p we consider the momenta {p^ i = 1, • • • , n} as auxiliarly variables, 
we can eliminate them on the basis of their equations of motion. One obtain 
the action 

So[q,u\ = j dtLo{q,q,u) (24) 

where the Lagrange multipliers u", a — 1, ■ ■ ■ ,m are seen now as real fields. 
The action ([24|) is invariant at the gauge transformations 

Sq^ = al,{q,q)s", Su" = s" - {V^" + GsUf,^)eP + {fl^ + Gsgl'^W . (25) 
The Noether's identities will have the form 

f^«^ + ^[-(^/f + TsU^) + + Gsgl-W] - |(gj) = 0. 

Starting from (IM)) we will develop the sp{2) BRST Lagrangean formalism [5]. 
The complete spectrum of the antifields is given in our case by 

<3^a = {<3Aa.<a} = {9*a>QLfc'^aa:.<a> ^,6=1,2}, (26) 
Qa = {QA^'^a} ^ {qi, Qaa,^a,Ua, a=l,2}. (27) 

It is well known that the Lagrangian dynamics is generated in a " anticanonical" 
structure. The generator of the Lagrangian BRST symmetry is 

S = S'o[g,w] H 
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Because So is unique one we will consider that S is unique and we will introduce 
two antibracket structures which have the same properties like in the standard 
theory: 

_ 6'F 5^G S^'F 6^G 

The functionals F and G are depend to Q-^ and Q\a- On the basis of the 
graduation properties and of the Grassmann parities [9] we define the paires 
canonical conjugate in respect with these antibrackets 

{Q*Ab,Q'')a^~5''^5,a (28) 

where are expressed by 1^5]) and are the fields of the theory (real q\ u" 
and ghosts Q"°,A") 

Q-^ EE {Q^, u"} {q\ A", u", a = 1, 2}. (29) 

For the Lagrange multiplyers we will have 

e(u") = £a,gh{u°') = 0. 

We observe that some antifields of the theory have a canonical conjugate in the 
antibracket structure (pS)) and other antifields not have canonical pair (P?]). So, 
the BRST differentials {s°, a = 1, 2} will have the following decomposition 

s% = (5°)'="' * = (*, Sf * +1/%, a = 1, 2 (30) 

where the non-canonical operators Va have the form [9] 

1/% = (_)-(Q^)e-''Q^,^ * . (31) 



'5Q, 



A 



The nilpotency condition for s° (|30p leads to the master equations ([15]). For our 
irreducible theory, the proper solution of the master eqs (|15p till terms linear in 
the antifields is 



+9,^(-)^n/^". + G5ff^^)(/;p + G5g^^)ecdg'"^g^^g'''^+ 
+uS ^uo. {-{Vp + + (-)^^(/|^ + G^^^)^^) A'^-f 
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+ Q.a i-r^' if p. + Gsgf,Z){X^Q^^ - + Gsg',;)e,,Q'^Q''^Q^'). (32) 

On the basis of the graduation rules and Grassniann parities [8], [9] we can 
identify the fohowing variables 

These identifications will be very useful in the gauge fixing procedure [TJ. On 
the basis of the identifications (155)) . we will have 

Y = r{q')u^. (34) 

The following relations are valid 

(5'' 1 S''f°'(Q') 
9« = J^i^^^^y^Y) = (35) 

'^^ = ^ ^ 

5"^ ,1 



. . {7^eabV,Y)=r{q^) (37) 

the remaining antifields vanishing because of the choice (p4|) for F. The gauge 
fixed action will be 

_ _ Sir fat i\ 

SiY ^ Si[Q ,u^^,Ua,u (q),g,^ = — u 



Sq'- 



aa 1 



_ 5'Y s[rm_ 

It leads to an effective action which is Sa-invariant and that can be further used 
in the path integral. This path integral can be written as 

^ J VQ'^Vul^VuaexpiiSiY)- (39) 
We will introduce the condensed notations 

cj,^ ^ {Q^, zi4 = {q\ g"^ r,K,,ua} (40) 

and we will consider the following BRST transformations 

where fj,^ are small fermionic constant parameters. The superjacobian of this 
transformations can be aproximate with supertrace (because it involves both 
fermionic and bosonic fields, the jacobian is a superdeterminant) and the new 
integrating measure will be 

Vcp^' =(^1- ((-1)-^^ + V^SW^ Vc^^ = (1 - iA^S)^,,W^. (42) 
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In the previous relation, the operators A°, a ~ 1,2 have the form (|16p . 
Exprimam A°-S using ([32]) : 



if^^ + Gs9';';)e,dQ'"'Q'''Q^')] - i-re^''ul,{^{^SabVbSY)KX'^ + 

-[f'^^^ + Gs9l^^)^{\eMY)}- 

The functions Gs which appear in the previous relations are the constraints 
in which the real momenta was substitute (on the basis of their equations of 
motions) in function oi q,q. In the same way appear the functions aj^ from 
{q\Ga.}. 

The antifields was eliminated by gauge fixing procedure (l55]) - ([57)) but con- 
sidering the gauge fixing functional of the form Y + SY. Then, A°-S = if 

and if -^{^eabVbSY) and g^{-^ScdVd6Y) anuleaza termenii din paranteza (can- 
cel the terms) ce apare pe lang With other words, the measure remain 
invariant (A^S* — 0) to the transformation (HTI) if 6Y induced by this trans- 
formation in the gauge fixing functional indeplineste cerintele de mai sus. If 
A^S ^ 0, the operators A", a = 1,2 determine the sp{2)- QME in the sp{2) 
symmetric formulation of gauge theories, where the action Sy+sy passes to a 
new "quantum action", W, containing the ?i— order corrections and satisfying 
the equations P^ . 



4 Conclusions 

The problem of the quantum extension of the BRST procedure has been inves- 
tigated, both in the BFV and in the BV formalisms. The main questions arising 
in our study were: What is the form of the quantum master equation? How is 
looking and what is the meaning of the quantum action? How is it possible to 
extend the standard BRST quantization procedure in order to obtain sp(2) or 
sp(3) quantum master equations? 
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We tried to offer comprehensible answers to all these questions and to il- 
lustrate them by considering the case of nonlinear gauge theories which mix 
bosonic and fermionic fields and whose gauge superalgebra is defined by non- 
linear relations of the form (|21I) . (P^ . For such a special system considered in 
|14] we implemented the Hamiltonian sp(2) formalism and we pointed out the 
condition for which the path integrals are BRST invariant and how the gauge 
fixing function has to be chosen in order to assure the invariance when it is not 
intrinsec for the model. 
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Abstract 

he paper takles with a procedure which allow to extend some linear, wave type models cur- 
rently used in describing phenomena appearing in atmosphere to the study of nonlinear models. 
More concretely, we present a practical way to generate the largest class of (1 + l)-dimensional 
second order partial differential equations (pdes) of a given form which could be reduced to an 
imposed ordinary wave type equation. This class generalize the ordinary differential equation 
describing the equatorial trapped waves generated in a continuously stratified ocean and will 
be obtained following the Lie symmetry and similarity reduction procedures. Moreover, some 
concrete nonlinear second order differential equations will be proposed as possible candidates 
for replacing more complicated, nonintegrable systems, as the Rossby type equation. 

Keywords; Nonlinear dynamical systems, Lie symmetries. Similarity reduction procedure, 
Rossby type symmetries. 

1 Introduction 

A rich variety of complex phenomena occuring in many physical fields, including the atmospheric 
dynamics, are described by linear differential equations which allow a simple handling of the con- 
straints which appear in the system's evolution. Although, the linearized models often do not 
adequately describe the dynamics of the processes as a whole and it is very simple to shift the 
system to a region in which the linear behavior is no longer valid. This is why, in order to capture 
the real behavior, to accurately estimate and control the complex systems in all their regimes, 
nonlinear models must be defined. In this case, the linear differential equations could appear as 
approximations to the nonlinear systems, valid under restricted conditions. 

The price to be paid when nonlinearity is taken into consideration appears in the investigation 
of the exact solutions of the attached equations. There are not standard methods of solving 
nonlinear differential equations, they are usually depend on the form of the equations and on their 
particular symmetries. Many interesting nonlinear models have been proposed over the last years 
[I] and a lot of methods have been developed in order to find solutions of equations describing 
these nonlinear phenomena. Some of the most important methods [7] are the inverse scattering 
method [8], the Darboux and Backlund transformations [9j, the Hirota bilinear method |10) . the Lie 
symmetry analysis [HI [111 ES] , etc. By applying these methods, many types of specific solutions 
have been obtained. For example, solitary waves or solitons, which have no analogue for linear 
partial differential equations, are very important for the nonlinear dynamical systems. 

In this paper we shall concentrate our attention to the Lie group method. It is well-known 
that this method is a powerful and direct approach to construct many types of exact solutions 
of nonlinear differential equations, such as soliton solutions, power series solutions, fundamental 
solutions [m [12], and so on. The existence of the operators associated with the Lie group of 
infinitesimal transformations allows the reduction of equations to simpler ones. The similarity 
reduction method for example is an important way of transforming a (1 + 1)— dimensional pde into 
an ordinary differential one. We shall concretely consider the inverse symmetry problem [16j and 
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we shall generate the largest class of second order (1 + 1)— dimensional pdes which generalize the 
ordinary, wave type, differential equation describing the equatorial trapped waves generated in a 
continuously stratified ocean. Practically, four types of waves appears in this case and have to 
be found among the solutions of the equation: Kelvin waves, Rossby waves, inertia-gravity waves 
and mixed Rossby-gravity waves. In the Boussinesq approximation and on an equatorial /3-plane, 
the equation which describe the m-th oscillation mode of the wave's vertical velocity (pmi^) the 
direction z has the form [23]: 

Q^2 + "Cl^'t'rniz) = (1) 

One consider for the velocity the boundary conditions: 

= —H) = (ocean floor) 
4>m{^ — 0) = (ocean surface) 

In the equation ([1]), Cm is a constant and Nm{z) represents the "buoyancy" frequency. Measure- 
ments made during El Nino events by 2 Japonese stations in the equatorial Pacific [23] show that 
the buoyancy frequency Nm{z) has strong variations with the water depth close to the surface and 
practically vanishes for higher depths. We notice that in the first case (at the surface, z G [0, 300] m) 
one can aproximate Nm{z) with an averaged value around N{z) = 2 • 10~^m • So, the equation 
([TJ can be linearized in one of the following forms: 

(f){z) = 0;z> 300 (2) 

-a N 

4>{z) + k^4>{z) = 0; z G [0, 300]; k = — = const. (3) 

In this paper, we shall consider the two previous wave type equations and we shall see how 
they can be extended towards (1 + 1)— second order differential equations with the same group of 
symmetry as the initial ordinary wave type equations have. 

The outline of this paper is as follow: after this introductory notes, in Section 2, we shall obtain 
the general determining system for a chosen class of (1 + 1)— dimensional models. The system will be 
generated by using the Lie symmetry approach and by asking for an imposed form of the similarity 
reduction equation. More exactly, we shall generate a class of (1 + 1)— second order differential 
equations which by similarity reduction come to the wave forms ([2]) and, respectively, Q. The 
general results of the second section will be particularized in Section 3, when concrete examples of 
two dimensional equations with similar solutions as the ordinary wave equations ([2]) and ([3]) will be 
generated. Moreover, we shall compute the form of the second order partial differential equation 
which admit an imposed form of symmetry, specific for the two dimensional Rossby type equation. 
So we shall be able to replace the study of this last strongly nonintegrable equation with a simpler 
class of equations observing similar symmetries. Some concluding remarks will end the paper. 



2 Determining equations for the Lie symmetry group 

Let us consider the class of general dynamical systems described in a (1 + l)-dimensional space 
(x,t) by a second order partial differential equation of the form: 

Ut = A{x, t)u2x + B{x, t)Ux + C(x, t)u Q{x, t, U, Ux,Ut,U2x) = (4) 

Our aim is to select from the general dynamical systems described by Q the class of differential 
equations which admit a similarity reduction to wave type equations of the form ([JD and As 
feed-back, the solution of the wave equations will be used in order to obtain a solution for Q. The 
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procedure that will be followed firstly implies to obtain the system of determining equations for 
the Lie symmetry group of (jj]). Then, an additional system of partial differential equations will 
be generated by imposing that ([!]) possess a reduced similarity equation of the wave type. Finally, 
this last system and the Lie determining equations will be solved and coefficient functions A{x,t), 
B{x,t), C{x,t) will be obtained. 

In this section we shall apply the Lie symmetry approach for the equation Let us consider 
a one-parameter Lie group of infinitesimal transformations: 

X = X + e^{t,x,u), t = t + eip{t, x,u), u = u + eri{t,x,u) (5) 

with a small parameter e ^ 1. The Lie symmetry operator associated with the above group of 
transformations can be written as follows: 

Odd 
U{x, t, u) = ip{x, t,u)—+ ^{x, + 7]{X, t, u)— (6) 

The second order equation Q{x,t,u, Ux,ut,U2x) = of the form ([H) is invariant under the action of 
the operator ([6]) if and only if the following condition [11] is verified: 

U^'H^) \n=o=0 (7) 

where U^"^^ is the second extension of the generator ([6]). A concrete computation shows that the 
coefficient functions from dH) and ([6]), A{x,t), B{x,t),C{x,t),if{x,t,u),^{x,t,u),r]{x,t,u), must 
satisfy the equation: 

if At + ^Ax)u2x + ivBt + iBx)ux + ipCtu + iCxU + Cr] + 5//^ - rf + At?^^ = (8) 

Coefficient functions r/^, rf^ rf^ appear in the process of extension of U towards f/*^^^ and their 
general expressions are given in |11) . Using these expressions in ([8]) and asking for the vanishing of 
the coefficients of each monomial in the derivatives of u{t,x), we obtain the following differential 
system: 

= 0; iPu = 0; = 0; V2u = 0; ^At + C^x + Aip^ - 2A^^ = 0; 

- ipBt - iBx + Bi., -it- B^t - 'iMxu + Ai2x = (9) 

-ifCtu - iCxU -Ct]- Br]^ + % + Crj^u - ipfiu - Ar]2x = 

The first four equations of the system Q lead, for coefficient functions (^(x, t, u),^ {x, t, u),r]{x, t, u), 
to the following reduced dependences: 

if = ip{t), e = C{x,t), 7] = M{x,t)u (10) 

Consequently, the remaining equations of dH) become: 

<fAt + ^Ax + A^t - 2Ai^ = 0; 

- ^Bt - CBx + BCx -it- B^t - '^AMx + AC2X = (11) 



-ipCt - iCx - BMx + Mt- iftC - AM2X = 

with 6 unknown functions: A{x,t),B{x,t),C{x,t) provided by the evolutionary equation and 
(p{t), ^(x,t), r]{x,t,u) introduced by the symmetry group of transformations ([5]) and described by 
the relations pop . 
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3 Similarity reduction procedure 

Let us consider now the similarity reduction procedure. In this section, some particular choices for 
the system ()lip will be considered. We shall find equations describing concrete dynamical systems 
which admit reduction through the similarity procedure to ordinary wave type equations of the form 
([2]) and ([3]). For the moment, we restrict the forms ([TO]) of the infinitesimals (p{t) , S,{x , t) , r](x , t, u) 
to the following separable expressions: 

if = ip{t), e = ^ix,t) = Ciix)^2it), V = M{x,t)u = Mi{x)M2{t)u (12) 

The Lie operator ([6]) becomes: 

U{x,t,u) = ^{t)-+ ^^{x)Ut)-Q^ + Mi{x)M2{t)u— (13) 

The general expressions of the invariants could be found if we should consider the characteristic 
equations associated with the new generator (|13p . These equations are: 



dt dx du 

W) " Ux)Ut) " M,{x)M2{t)u ^ ^ 

By integrating the previous equations, two invariants are obtained with the following expressions: 

In the similarity reduction procedure two similarity variables have to be considered: 

h = z, l2 = Hz) (16) 

The invariants (fT6]) allow us, by an appropriate change of coordinates {u,x,t} — )• {z,(l){z)}, to 
reduce the initial (1 + 1) dimensional equation to an ordinary differential equation of the form: 

n'[z,Hz),Hz),-] = o (17) 

3.1 Homogeneous wave type equation 

Our aim is now to select from the general dynamical systems described by dH the class of differential 
equations for which the equation ()17p can be reduced at a wave type equation of the form ([2]): 

= 0^^(z) = az + 6 (18) 

where a and b are arbitrary constants. 

The previous solution, written in terms of the initial variable (x,t), leads to the following form 
of the solution it(a;, t) of ([!]): 



u(x, t) 



aexp ( / dx — f ^ } dt ] + b 

Ci(x) J (fit) 



(19) 



For convenience reasons, we shall impose the following relations to be valid: 
with g, V arbitrary constants. 
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In terms of notations ()20p . the infinitesimals ()12p and the solution (|19p become: 



u{x, t) 



7? 



qv ^— — u 



(21) 
(22) 



■P(x)' ' ^ P(a;) 
[a exp {P{x) — qt) + 6] exp {vR{x)) 

The solution (|22p must verify the equation ^ which describes the analyzed model. This condition 
generates a differential system of the form: 

= q + 2vA{x,t)R{x)P{x) + v'^A{x,t)[R{x)f + A{x,t)P{x) + A{x,t)[P{x)f + 

+vA{x, t)R{x) + B{x, t)P{x) + vB{x, t)R{x) + C{x.t) 
= v'^A{x,t)[R{x)f + vA{x,t)R{x) + vB{x,t)R{x) + C{x,t) (23) 

For an unitary analysis, it is necessary to describe the differential system (|lip . obtained in 
the previous subsection, in terms of the functions P{x) and R{x) introduced by (pO|) . Using the 
expressions (j2ip we obtain the following differential system: 



ipAtP'^ + qipA^P^ + iptAP^ + 2qipAP2x 







ipBtP^ + q^B^P^ + qipBP2.P^ + qip^P! + iPtBP^+ 
+ 2vq^AR2^Pl - 2vqipAR^P2^P^ + q^AP^^P^ 



(24) 



2qipAP,Pl = 

^CtP^ + q^C.P^ + qv^BR2^Pl - qvipBR,P2,P^ + ^^CP^ + vq^ARs.P^- 

+vqipAR^P3^P^ - 2vq^AR2APl + +2qvipAR^P^Pl - qv^^R^P^ = 

Conclusion: Our problem is to find the class of (1 + 1) evolutionary equations of type (HD 
which could be reduced by the similarity approach to an ordinary wave type equation. Solving this 
problem is equivalent with searching the solutions of the system described by equations (|23p and 

(EH). 

Remark 1: The system ()23p -()24 p can be solved following two paths: (z) by choosing a concrete 
dynamical system, that is to say concrete expressions for the functions A{x,t)^ B{x,t), C{x,t) and 
trying to find out if this equation admits or not solution of the type (|22p . Now the unknown 
functions of the system are ip{t), P{x), R{x) defined by (pO|) : {ii) by considering A{x,t), B{x,t), 
C{x,t) as unknown functions and by choosing (p{t), P{x), R{x). This is the way we shall follow in 
the next section. 

Remark 2: In the case (ii) the general solutions obtained by computational way can be 
expressed as: 

'qt - P{x)' 



A{x,t) 
B{x,t) 
C{x,t) 



-2F 



qt-P{x) 



exp [G(x)] 

i[P(x)]2 + vR{x)P{x) + ^P{x)) exp [-G{x) 



qt-P{x) 
1 



P{x) 

R{x)P{x) + P{x){-R{x) + R{x){vR{x) + P{x))) \ exp[-^(^)l +qR{x) 



where 



G{x) 



P{x) 



2{D(^y){P){a)^ ( PM±^^P(-1) q + [D{P){a)f D{ip) ( ^W+g^-^W 



DiP)ia)^ 



P{a)+qt-P(x) 



da 



(25) 



(26) 



These solutions are valid for arbitrary constants q, v and for an arbitrary function F 



qt-P(x) 
1 
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3.2 Harmonic oscillators 



Let us consider now that the similarity reduction equation is an ordinary oscilator type equation 
of the form 

'^ + k^iz) = (27) 

It has the solution: 

(j){z) = asm{kz) + bcos{kz) (28) 

Here k, a and b are arbitrary constants. 

To writte down the previous solution in terms of the initial variable u{x, t) means that the 
solution of dl]) should have the form: 

u(x,t) = aexp I f ^^^^} dx] sin (Vkexpl f -^dx - f ^^^dt]] + (29) 

J Ciix) J \ \J ^i(x) J ip{t) J J 

+6exp ( —M-r- f , } dx] cos ( Vkexp ( f dx — f ^^ j' } dt 



For convenience reasons, we shall impose again the following relations to be valid: 



-Pit) m J (i(x) " ' "J {lOc 

with q, V arbitrary constants. 

In terms of notations (|3Up . the infinitesimals (|12p and the solution (j29|) become: 

^ = ^[t), C = Q-^, V = A, X ^ (31) 
F[x) F[x) 



u{x,t) = [asm{kexp {P{x) — qt)) + bcos{kexp {P{x) — qt))]exp {vR{x)) 

(32) 

The solution ()32p must verify the equation ^ which describes the analyzed model. This condition 
generates the vanishing of the coefficient function A{x,t) and two other differential equations of 
the form: 

A{x,t) = 

q + B{x,t)P{x) = (33) 
vB{x,t)R{x)+C{x,t) = 

For an unitary analysis, it is again necessary to describe the general differential system (|lip 
obtained in the previous section, in terms of the functions P{x) and R{x) introduced by ([30]) . 
Taking into account the equations (|3ip . we obtain the following differential system: 

= ipBtP^ + qipB,P^+qipBP2.P^+q^tP! + iptBP^ 
= ifCtP^ + qipC^Pl + qvipBR2^Pl - qvipBR^P2,.P^ + 

+VtCP^ - qv^tRxPl (34) 

The system (|33p -(|34p can be solved following two paths: (z) by choosing a concrete dynamical 
system, that is to say concrete expressions for the functions B{x,t), C{x,t) and trying to find out 
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if this equation admits or not solution of the type ()32p . Now the unknown functions of the system 
are f{t), P{x), R{x) defined by ([30]) : (ii) by considering B{x,t), C{x,t) as unknown functions and 
by choosing ip{t), P{x), R{x). 

This second case is the way we are interested in to fohow and, in this case, the general solutions 
obtained by computational way can be expressed as: 

B(..t) = ^.C(.,*) = ^ (35) 

or in terms of the coefficient functions v'(t), £,{x^t), M{x,t) which appear in the general Lie sym- 
metry operator (??), in the equivalent forms: 



fit) ' ^ ' ^ V{t) 
3.3 Rossby type symmetries 

The equation for coupled gravity, inertial and Rossby waves in a rotating, stratified atmosphere 
using the /3-plane approximation (which simplifies the spherical geometry whilst retaining the 
essential dynamics) and the Boussinesq approximation which filters out higher frequency acoustic 
waves can be written in (2 + 1)— dimensions in the form 



dt V9x2 + Qy2 )Uy--P (36) 

As we mentioned, this equation describes the coupling between the inertial, the gravity and the 
Rossby waves, but also the shallow water in an ocean of depth H. It was proven |25j that reducing 
the model to (1 + 1)— dimensions, (x,t), it admits some very simple Lie symmetries of the form: 

ip{t) = ct + ci, ^(x, t) = cx + f{t), rj{u) = —3cu (37) 

with f{t) arbitrary function and cq, c arbitrary constants. Despite this fact, it is still difficult to 
find explicite solutions for the equation p6p . This is why, we shall consider another approach: we 
shall impose the Lie symmetries ()37p to our general equation @ and we shall try to find the class of 
equations which observe them. This means that we have in fact to impose the Rossby symmetries 
P7p to the system (|lip . It will take the form: 



{ct + ci)At + {cx + C2)A^ + 3cA = 
-{ct + ci)Bt- {cx + C2)Ax-2cB = 

-{Ct + Cl)Ct - {cx + C2)Cx - cC = 



with the unknown functions A{x,t), B{x,t), C{x,t). 
This system admits the solutions: 



A{x,t) 
B{x,t) 
C{x,t) 



F{x{ct + ci) - 


C2t) 


F{xip{t)- 


C2t) 


{Ct + Ci)3 








G {x{ct + ci) - 


C2t) 


G {xip{t) - 


C2t) 


(Ct + Ci)2 








H{x{ct + ci) - 


-C2t) 


H{x(f{t) - 


-C2t) 


ct + Cl 




ifit) 





with F, G, H arbitrary functions of their arguments. These expressions give us equations of the 
form (jH which are equivalent from the point of view of their symmetries with the Rossby equation. 
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4 Conclusions 



The problem of finding exact solutions for nonlinear differential equations plays an important role in 
the study of nonlinear dynamics. There are many ways of tackling with it. One of them is based on 
the Lie symmetry method. This method supposes to find the symmetries of the system and, on this 
basis, to try to determine the general or some particular solutions of the equations. There is a direct 
approach in which the symmetries of a given equation are obtained, but also an inverse problem has 
been formulated [16] . A step forward for this latter approach is represented by the use of similarity 
reduction, a procedure which allows the reduction of the number of degrees of freedom and, by 
that, simplifies the problem of solving the equation. This paper used this approach and determined 
a class of (1 + 1) dimensional second order diff'erential equations which can be reduced to ordinary 
wave-type equations with simple solutions. Using the Lie symmetry and the similarity reduction 
procedures, some particular cases of the equation ^ arise as good candidates of equations which 
could be used as generalization of the linear wave type equations describing complex atmospheric 
phenomena. Moreover, following our method, we were able to write down the solutions of these 
equations, solutions which otherwise could be derived by computational methods, but in a very 
complicated form. Another interesting results of our paper consisted in the fact that a complicated, 
nonintegrable equation, the Rossby equation, could be replaced by another, simpler equation, which 
have similar symmetries. The paper is important both by these results, but also as a methodological 
approach in finding exact solutions through similarity reduction procedure. We have shown how, 
starting from a particular form of solution for the reduced equation we could recover the solution 
of a most complicated problem, defined in a space with more than one dimensions. We tackled 
out a particular case, looking only for linear solutions of the reduced equation and considering that 
the coefficient functions appearing in the symmetry operators are separable. The problem can be 
extended for other cases, too. 

Acknowledgements: This work was supported by GNOSIS -UEFISOSU, project number PNII 
- IDEI code 418/2008. 
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